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ABSTRACT 

Galaxy clusters possess turbulent magnetic fields with a dominant scale length /g ~ 
1 — 10 kpc. In the static-magnetic-field approximation, the thermal conductivity Kj for 
heat transport over distances 3> /g in clusters is ~ Ks/g/Ls(Pe), where Ks is the Spitzer 
thermal conductivity for a non-magnetized plasma, the length Ls (ro) is a characteristic 
distance that a pair of field lines separated by a distance ro < /g at one location must 
be followed before they separate by a distance /g, and p e is the electron gyroradius. 
We introduce an analytic Fokker-Planck model and a numerical Monte Carlo model 
of field-line separation in strong magnetohydrodynamic (MHD) turbulence to calcu- 
late Ls(ro). We also determine Ls( r o) using direct numerical simulations of MHD 
turbulence with zero mean magnetic field. All three approaches, like earlier models, 
predict that Ls asymptotes to a value of order several /g as ro is decreased towards Id 
in the large-Zg/Z^ limit, where Id is the dissipation scale, which is taken to be the 
proton gyroradius. When the turbulence parameters used in the Fokker-Planck and 
Monte Carlo models are evaluated using direct numerical simulations, the Fokker- 
Planck model yields Ls(p e ) — 4.51b and the Monte Carlo model yields Ls(p<?) — 6.5/g 
in the large-Zg/Z^ limit. Extrapolating from our direct numerical simulations to the 
large-/g//^ limit, we find that Ls(p e ) — 5 — 10/g, implying that Kj — O.lKs — 0.2ks 
in galaxy clusters in the static-field approximation. We also discuss the phenomenol- 
ogy of thermal conduction and particle diffusion in the presence of time-varying tur- 
bulent magnetic fields. Under the questionable assumption that turbulent resistivity 
completely reconnects field lines on the time scale lg/u, where u is the rms turbulent 
velocity, we find that Kj is enhanced by a moderate amount relative to the static-field 
estimate for typical cluster conditions. 
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1. Introduction 



In the cooling-flow (CF) model of intracluster plasma, radiative cooling causes plasma to 
flow towards a cluster's center and cool to sub-x-ray temperatures, presumably ending up as ei- 
ther stars, smaller compact objects, and/or cold gas. Aside from the gravitational work done on 
inflowing plasma, heating of intracluster plasma is neglected in the model, and mass accretion 
rates M are as high as 1O 3 M yr _1 for some clusters (Fabian 1994). A longstanding problem for 
the CF model has been the difficulty in accounting for all the accreted mass. For example, the 
observed rates of massive star formation are a factor of 10-100 less than expected if the cooling 
plasma predicted by the model ends up forming stars with a normal IMF (Crawford et al 1999, 
Fabian 2002). In addition, recent x-ray observations find no evidence of plasma cooling to tem- 
peratures below 1-2 keV (Peterson et al 2001, Tamura et al 2001). These difficulties suggest that 
some form of heating approximately balances radiative cooling, thereby dramatically reducing M 
relative to CF estimates. A number of heating mechanisms have been considered, such as galaxy 
motions (Bregman & David 1989), supernovae (Bregman & David 1989), cosmic-rays (Bohringer 
& Morfill 1988, Tucker & Rosner 1983), active galactic nuclei (Pedlar et al 1990, Tabor & Binney 
1993, Binney & Tabor 1995, Ciotti & Ostriker 2001, Churazov et al 2002), dissipation of turbulent 
energy (Loewenstein & Fabian 1990, Churazov et al 2003, Chandran 2003), and thermal conduc- 
tion, which can transport heat from the hot outer regions of a cluster into the relatively cooler core 
(Binney & Cowie 1981, Tribble 1989, Tao 1995, Chandran & Cowley 1998, Narayan & Medvedev 

2001, Gruzinov 2002, Voigt et al 2002, Zakamska & Narayan 2002). For thermal conduction to 
approximately balance cooling, the thermal conductivity must be a significant fraction of the 
Spitzer value for a non-magnetized plasma, Ks, and in some clusters even greater than Ks (Fabian 

2002, Zakamska & Narayan 2002), where 



Coulomb logarithm (Spitzer 1962). 

Galaxy clusters are filled with tangled magnetic fields with a dominant length scale Is — 
1 — 10 kpc (Kronberg 1994, Taylor et al 2001, 2002) that is much less than the size of a cluster 
core, R c ~ 100 kpc. Both optical-line-emitting gas in clusters and hot intracluster plasma are 
observed to be in turbulent motion (Fabian 1994, Churazov et al 2003). The effects of turbulent 
magnetic fields and velocities on K T are the subject of this paper, and have been investigated by a 
number of authors (e.g., Tribble 1989, Tao 1995, Chandran & Cowley 1998, Chandran et al 1999, 
Narayan & Medvedev 2001, Malyshkin & Kulsrud 2001, Gruzinov 2002). 



In terms of this definition, in which Ks is expressed as a diffusion coefficient, the heat flux is given by — n e fcflKs VT. 




(1) 



T is the temperature, kg is the Boltzmann constant, n e is the electron density, and In A is the 
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Transport in the presence of strong turbulence is a difficult and unsolved problem. It is likely 
that the thermal conductivity for a particle species scales like the test-particle diffusion coefficient 
for that particle species (Rechester & Rosenbluth 1978, Krommes, Oberman, & Kleva 1983): 
since the collisional transfer of energy between particles occurs locally in space, diffusion of heat 
accompanies the diffusion of heat-carrying particles. We thus estimate Ky in clusters from the 
relation 

— — — , (2) 
K S D ' 

where D is the diffusion coefficient of thermal electrons and Do is the thermal-electron diffusion 
coefficient in a non-magnetized plasma. 234 

Since thermal electrons in clusters move much faster than the E x B velocity of field lines, 
a reasonable first approximation for D is obtained by treating the magnetic field as static. In a 
collisional plasma, particle diffusion over distances ^> l B in a static field depends critically on the 
rate of separation of neighboring magnetic field lines (Rechester & Rosenbluth 1978, Chandran 
& Cowley 1998). If two field lines within a snapshot of strong magnetohydrodynamic (MHD) 
turbulence are at some location separated by a distance tq <C l B , they will separate by a distance l B 
after some distance z along the magnetic field. The length Ls (ro) is a characteristic value of z 
defined in section 2. In the static-magnetic field approximation, the thermal conductivity K T in 
galaxy clusters over distances l B satisfies 

Ls(pe) 

where p e is the electron gyroradius. Equation (3) makes use of the fact that electron motion along 
the magnetic field is relatively unimpeded by magnetic mirrors and whistler waves excited by the 
heat flux because the Coulomb mean free path is short compared to both l B and the temperature- 
gradient length scale (see section 3). The factor l B /Ls(p e ) in equation (3) measures the reduction 
in Kt associated with tangled field lines, which increase the distance electrons must travel in going 
from hotter regions to colder regions. 



2 Electrons generally make the dominant contribution to Kj because they diffuse more rapidly than ions. 

3 In a steady state, an electric field is set up to maintain quasi-neutrality. This electric field reduces by a factor 
of ~ 0.4 in a non-magnetized plasma (Spitzer 1962). We do not consider the effects of turbulence on this reduction 
factor. 

4 There is some ambiguity in the right-hand side of equation (2). We take D and Do to be diffusion coefficients for 
particles of a specified energy as opposed to diffusion coefficients of a particle that diffuses in energy as well as space. 
Thus, D/Do may scale differently for electrons of energies, e.g., kgT and IkgT. In the static-field approximation this 
is not an issue since D/Dqis the same for all energies of interest for cluster parameters. However, when we estimate 
the effects of turbulent resistivity, we pick a representative energy at which to evaluate D/Dq as discussed in section 6. 
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Three previous studies (Jokipii 1973, Skilling, Mclvor, & Holmes 1974, Narayan & Medvedev 2001) 
calculated Ls(ro) for strong turbulence assuming a power spectrum of magnetic fluctuations on 
scales ranging from /g to a much smaller dissipation scale lj, which is approximately the proton 
gyroradius p ; in galaxy clusters (Quataert 1998). Using different definitions of Ls, each study 
found that Ls(ro) is of order Zg in the small-(ro//g) and small-Lj/Zg limits provided ro is not vastly 
smaller than lj. Jokipii (1973) assumed an isotropic turbulent magnetic field and employed a 
stochastic model of field-line separation, which he solved using a Monte Carlo numerical method. 
Skilling et al (1974) assumed isotropic turbulence and introduced an approximate equation of the 
form d(r)/dz = F((r)) for the average separation r of a pair of field lines, with the function F 
estimated from the power spectrum of the turbulence. Narayan & Medvedev (2001) introduced 
a similar equation for the evolution of the mean square separation, d(r 2 ) /dl = G((r 2 )), and esti- 
mated G using the Goldreich-Sridhar (1995) model of anisotropic MHD turbulence. 

In this paper, we calculate Ls using three different methods. First, we consider an approximate 
model in which the separation r of a pair of field lines evolves stochastically and is described by 
the Fokker-Planck equation 



where y = ln(r/ /g), / is distance along the field in units of /g, Ay is the increment to y after a field- 
line pair is followed a distance Al along the field, (...) is an average over a large number of field- 
line pairs, and P(y\,l)dy is the probability that y is in the interval (yi,y\ +dy) after a distance / 
along the field. Our model is similar to Jokipii's (1973), although we determine the functional 
form of (Ay)/ A/ and ((Ay) 2 /A/ using the Goldreich-Sridhar model of locally anisotropic MHD 
turbulence, and we solve the Fokker Planck equation analytically, which allows us to determine the 
functional dependence of Ls on turbulence parameters. If we evaluate the turbulence parameters 
in the model using data from direct numerical simulations, we find that Ls(ro) — > 4.5/g as ro is 
decreased towards l d in the large-/g//^ limit, and that Ls(pe) — ^s(^)- 

Our second method for calculating Ls (ro) uses a numerical Monte Carlo model of field line 
separation in which each random step Ay is of order unity (the Fokker-Planck equation assumes 
infinitesimal Markovian steps). When the model parameters are evaluated using data from direct 
numerical simulations, the Monte Carlo model gives Ls(p e ) ~ 6.5/g in the large-/g/Lj limit. 

Our third method for calculating Ls(ro) involves tracking field-line trajectories in direct nu- 
merical simulations of MHD turbulence with zero mean magnetic field. The results of our highest 
resolution simulations are consistent with the prediction of the theoretical studies that Ls (ro) ap- 
proaches a value of order several /g as ro is decreased towards l c i in the large-/g / lj limit. Extrapo- 
lating our numerical results to the large-Zg/Zj limit suggests that Ls(p e ) ~ 5 — 10/g in clusters. 

Field evolution and turbulent fluid motions increase electron (and ion) mobility, enhancing K T 




(4) 
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to some degree. Turbulent diffusion in clusters has been studied by Cho et al (2003). In this paper 
we develop a phenomenology to describe the interplay between field evolution and single-electron 
motion under the questionable assumption that turbulent resistivity completely reconnects field 
lines on the time scale Ib/u, where u is the rms turbulent velocity. A similar assumption was 
explored by Gruzinov (2002). We find three limiting cases for the thermal conductivity. For 
typical cluster parameters and turbulent velocities (Churazov et al 2003) Kj ~ \/k^uIb, a value 
that is somewhat larger than both the turbulent diffusivity ~ ul B for that cluster and the static- 
field estimate of Ks/g/Ls ~ O.Iks — 0.2ks. More work, however, is needed to clarify the role of 
turbulent resistivity before firm conclusions can be drawn about its effects on Kj. Additional work 
is also needed to quantify factors of order unity that have been neglected in estimating both in 
the presence of efficient turbulent resistivity and in the static-field approximation. 

The remainder of this paper is organized as follows. In section 2 we review the phenomenol- 
ogy of thermal conduction in static tangled magnetic fields. In section 3, we discuss the effects 
of magnetic mirrors and microturbulence on electron diffusion along field lines. We present our 
Fokker-Planck and Monte Carlo models of field-line separation in section 4. We compare these 
approximate theoretical models with results from direct numerical simulations in section 5. We 
estimate the effects of turbulent resistivity on in section 6 and summarize our conclusions in 
section 7. We present numerical simulations of field-line separation for different types of MHD 
turbulence in a companion paper (Maron, Chandran, & Blackman 2003). Table 1 defines some 
frequently used notation. 



2. The phenomenology of thermal conduction in a static tangled magnetic field 

We assume that the magnetic fluctuations possess an inertial range extending from an outer 
scale Ib to a much smaller inner scale 1^ with the magnetic energy dominated by scales ~ /g. 
Except where specified, the discussion focuses on the case relevant for clusters in which the mean 
magnetic field is negligible. 

A tangled magnetic field line is essentially a random-walk path through space. If a particle is 
tied to a single field line and travels a distance / ^> l B along the static magnetic field, it takes ~ I/Ib 
random-walk steps of length ~ Is, resulting in a mean-square three-dimensional displacement of 

((A;c) 2 ) = a/ fi /, (5) 

where a is a constant of order unity (values of a for the numerical simulations used in this paper 
are listed in table 2). When there is a mean field Bo comparable to the rms field, ((A*) 2 ) in 
equation (5) is interpreted as the mean-square displacement perpendicular to Bq. If the particle's 
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Notation Meaning 



l B dominant length scale of the magnetic field 

l d magnetic dissipation scale 

(/) average distance in units of /g that a field-line pair must be followed before separating by a 
(I 2 ) average square of the distance in units of l\ that a field-line pair must be followed before separatir 
(ro) h {I 2 ) I (I) — characteristic distance a pair of field lines separated by a distance tq must be followed b( 

p e electron gyroradius 

Kj thermal conductivity 

Ks Spitzer thermal conductivity in a non-magnetized plasma 

D 3D single-electron diffusion coefficient 

Z>0 3D single-electron diffusion coefficient in a non-magnetized plasma 

D|| diffusion coefficient for electron motion along the magnetic field 

X thermal-electron Coulomb mean free path 

u rms turbulent velocity 



Table 1: Definitions. 

motion along the field is diffusive with diffusion coefficient Dy, then 

l ~ y^y , (6) 

and (Rechester & Rosenbluth 1978, Krommes et al 1983) 

<(A*) 2 )oc' 1/2 , (7) 
indicating subdiffusion: D = lim((Ax) 2 )/6? — > (see also Qin et al 2002a,b). This process is 

t — >oo 

called double diffusion: the particle diffuses along the field line, and the field line itself is a random- 
walk path through space. 

The vanishing of D for a particle tied to a single field line can be understood by considering 
a particle starting out at point P in figure 1 . If this particle moves one Coulomb mean free path X 
along its field line Fi towards point Q and then randomizes its velocity due to collisions, it has 
a ~ 50% chance of changing its direction of motion along the magnetic field and returning to 
its initial point R In contrast, in a Markovian three dimensional random walk, the second step 
is uncorrelated from the first and there is a vanishing probability that a particle will return to its 
initial location. In a cluster, pe/h — 10~ 15 , and thus it is tempting to assume that electrons are 
tied to field lines. The importance of equation (7) is that any study that assumes that electrons are 
perfectly tied to field lines will conclude that K T = D = 0. 
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Of course, an electron is not tied to a single field line. As pointed out by Rechester & Rosen- 
bluth (1978), small cross-field motions enhanced by the divergence of neighboring field lines leads 
to a non-zero D. This can be seen with the aid of figure 1. Suppose an electron starts out at 




Fig. 1. — Trajectory of a diffusing electron. 

point P on field line Fi traveling towards point Q. After moving a short distance, field gradients 
and collisions cause the particle to take a step of length ~ p e across the magnetic field, from Fi 
to a new field line F2. Although the electron continuously drifts across the field to new field lines, 
let us assume for the moment that it remains attached to F2. As the electron follows F2, F2 di- 
verges from Fi . Let z be the distance that F2 must be followed before F2 separates from Fi by a 
distance (Because the electron continuously drifts across the field, it typically separates from 
Fi after traveling a distance somewhat less than z along the field; this effect, however, is ignored in 
this paper.) After the electron moves a distance z along F2, its subsequent motion is not correlated 
with F\ . The electron proceeds to point R, and then its collisional random walk along the magnetic 
field changes direction, bringing it back towards point Q. Instead of following F2 back to point Q, 
however, the electron drifts across the field and ends up on a new field line F3. After following F3 
for a distance ~ z, the electron separates from F2 by a distance ~ Is and proceeds to point S. 

In this example, the electron's small cross-field motions and the divergence of neighboring 
field lines allow the electron to escape from its initial field line and undergo a Markovian random 
walk in three dimensions. The fundamental random walk step is a displacement of length mz along 
the magnetic field, where m is some constant of order unity, perhaps 2 or 3. From equation (5), 
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when mz ^> Ib, a single random step corresponds to a 3D displacement of 

(A*) 2 ~ amzl B - (8) 

[When mz ^> Ib, the difference between the actual value of (Ax) 2 and its mean becomes small.] 
When mz 3> X, where X is the Coulomb mean free path, a single step takes a time 

2 2 

m z 

Af~-^-. (9) 

When mz is only moderately greater than Ib or X, equations (8) and (9) remain approximately valid. 
During successive random steps, a particle will find itself in regions of differing magnetic shear, 
and thus z will vary. The diffusion coefficient is given by D = ((Ax) 2 )/6{At) where (. . .) is an 
average over a large number of steps (Chandrasekhar 1943). Ignoring factors of order unity we 
obtain 5 

D\\l B 

D~-^, (10) 
Ls 

with 

(z 2 ) 
(z) 

In general, Ls, (z), and (z 2 ) are functions of the initial separation of a field-line pair, tq, but for 
electron thermal conduction we have set tq = p e . If there is a mean magnetic field comparable 
to the fluctuating field, equation (10) is recovered provided D is replaced by D^, the coefficient 
of diffusion perpendicular to the mean field. We note that taking Ls = (z) leads to similar results 
since (z) ~ (z 2 ) / (z) in our direct numerical simulations and Fokker-Planck and Monte Carlo calcu- 
lations. In clusters Dm ~ Do as discussed in section 3, where Do is the electron diffusion coefficient 
in a non-magnetized plasma. As a result, equations (2) and (11) give 

*r - (12) 
Ls 



3. Electron diffusion along the magnetic field 

For mono-energetic electrons subject to a fixed Coulomb pitch-angle scattering frequency, 
the diffusion coefficient Dy for motion along a uniform magnetic field is equal to the three- 
dimensional diffusion coefficient Do for motion in a non-magnetized plasma. Two mechanisms 



5 In applying equation (8) first and then averaging z and z 2 , we are averaging separately over the wandering of a 
single field line through space and the separation of neighboring field lines. This is justified to some extent since the 
former is dominated by the value of B at the outer scale and the latter depends on the magnetic shear throughout the 
inertial range. Although some error is introduced by averaging separately, equations (10) and (11) are sufficient for 
estimating K^. 
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for suppressing relative to have been discussed in the literature: magnetic mirrors and 
wave pitch angle scattering. Magnetic mirrors associated with a cluster's turbulent magnetic field 
significantly reduce D\\ only when X > /g, where X is the Coulomb mean free path of thermal 
electrons (Chandran & Cowley 1998, Malyshkin & Kulsrud 2001). In cluster cores, however, X 
is significantly less than /g, and thus mirrors have only a small effect. [For Hydra A and 3C295, 
X ~ 1 kpc at 100 kpc from cluster center, and X — 0.02 — 0.05 kpc at 10 kpc from cluster center 
(Narayan & Medvedev 2001).] When the Knudsen number Nk = X/L T \\ approaches 1, where 
L r || = T/\b ■ VT\ and b is a unit vector pointing along the magnetic field, the heat flux becomes 
large and excites whistler waves that enhance the pitch-angle scattering of electrons and reduce 
(Pistinner & Eichler 1998) 6 . However, for heat conduction into a cluster core, |Vr| ~ T/R c with 
R c ~ 100 kpc, and Nk "C 1- Thus, wave pitch-angle scattering has only a small effect, and 

D||~D . (13) 

We note that the chaotic trajectories of field lines cause L T ^ to be larger than T/|Vr|. 



4. Fokker- Planck and Monte Carlo models of field-line separation in strong MHD 

turbulence 

We adopt the Goldreich & Sridhar (1995) model of strong locally anisotropic MHD turbu- 
lence, which is supported by direct numerical simulations (Cho & Vishniac 2000, Maron & Gol- 
dreich 2001, Cho & Lazarian 2003). We assume that the fluctuating field is equal to or greater than 
any mean field in the system. The separation of neighboring magnetic field lines in strong MHD 
turbulence is dominated by shear Alfven modes. On scales smaller than /g, an Alfven-mode eddy 
is elongated along the direction of the average of the magnetic field within the volume of the eddy, 
denoted Bi oca i, with (Goldreich & Sridhar 1995, Cho & Vishniac 2000, Maron & Goldreich 2001, 
Lithwick & Goldreich 2001) 

fX \ 1/3 

5x ± ~5iocaif-^J , (14) 

and 

X^xfl l J\ (15) 

where is the rms magnetic fluctuation of an Alfven-mode eddy of width X^ measured across 
Biocai and length A,|| measured along Bi oca i. In fully ionized plasmas, the dissipation scale Id 
for Alfven modes is set by collisionless effects, and is comparable to the proton gyroradius p ; 



6 The formula in Pistinner & Eichler's (1998) paper is k t /k s - 1/[1 +250p e (A,/L r .|| )], but the factor of 250 should 
be corrected to a factor of 10 (Pistinner, private communication) 
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(Quataert 1998). The magnetic-field perturbation of an Alfven mode is perpendicular to Bi oca i. 
Equations (14) and (15) thus imply that when two field lines separated by a distance r are fol- 
lowed for a distance r 2 / 3 /^ 3 , r either increases or decreases by a factor of order unity assuming 
Id < r < h (Narayan & Medvedev 2001, Maron & Goldreich 2001, Lithwick & Goldreich 2001). If 
r < Lj, the separation or convergence of the field lines is dominated by the eddies of width lj, and r 

2/3 1/3 

increases or decreases by a factor of order unity when the field lines are followed a distance I J l B ' 
(Narayan & Medvedev 2001). We define 

/ (r// s ) 2 / 3 m d <r<l B 
Al= {(l d /l B )^ ifr<l d ' (16) 

y = ln(r/l B ), (17) 
Ay to be the change in y when the field lines are followed a distance Al, 

a=(Ay), (18) 



with a taken to be positive, and 



where a and b are of order unity. 



b=((Ay) 2 )/2, (19) 



To obtain Monte Carlo and analytic solutions for Ls, we make several approximations. The 
changes in r over a distance r 2 / 3 /^ 3 along the field associated with eddies of width much smaller or 
larger than r (or l c { if r < l c {) are small compared to the changes arising from eddies of width r (or l d 
if r < Id) and are neglected. We also take a and b to be independent of / and y, and consecutive 
values of Ay to be uncorrected. To obtain an approximate analytic solution for Ls, we make the 
further approximation of describing the stochastic variation of y with the Fokker-Planck equation 



dP_ d_ 
~dl~~dy~ 



Al 



P 



32_ 
+ dy 2 



((Ay) 2 ) 
2AZ 



P 



(20) 



where P(y',l)dy is the probability that y is in the interval (/,/ + dy), and / is distance along 
the magnetic field in units of Zg. The additional approximation in introducing equation (20) is 
associated with y changing by order unity during a single random step [equation (20) assumes 
infinitesimal steps]. 

We now solve equation (20) to obtain an analytic solution for Ls. Substituting equations (16), 
(18), and (19) into equation (20) yields 
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where 



ae -2y/3 p _ |_ { be -2yP P \ if yd < y < o 

r=< ^ V , (22) 

ae-^p _ be - 2 yj/i^- if y<y d 



y,i 1" ( -r- ) • (23) 



is the probability flux, and 

v.j = In I 

which is the value of y at the dissipation scale. We solve equation (21) with initial condition 
P(y) = S(y — yo) at / = and boundary conditions i> = 0aty = 0, P^Oasy^ and P and 
r continuous at y = yd- For electron thermal conduction in galaxy clusters, the quantity of interest 
is Ls when the initial separation r$ is the electron gyroradius, and thus we take yo < yd- The 
boundary condition P = at y = means that T(l,y = 0)dl gives the probability that a field-line 
pair separates to a distance lg for the first time after a distance between / and / + dl along the field. 

We proceed by making the substitution 

P = x m f, (24) 

with 

x = e y l 3 (25) 

and 

m = 2 + 3a/2b. (26) 
We then take the Laplace transform of equation (21), with the Laplace transform of / defined by 

/'CO 

7(5) = / f(l)e~ sl dl. (27) 
Jo 

For %d < x < 1, where Xd = e yd ^ 3 is the value of x at the dissipation scale, 

a 2 7 137 /v 2 9s\- n 

J +-J--l- + -)f = 0, (28) 



dx 2 x dx \ x 2 



with 



Since /(l) =0, 



3a 



(29) 



/ = d [/ v (\|/x)# v (\|/) -# v (\|/x)/ v (\|/)], (30) 
where 7 V and £y are modified Bessel's functions, 
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and c\ is a constant to be determined by applying the boundary conditions at x d after the solution 
for / for x < x c i has been obtained. 



For x < x d , 



3 2 / , 5 8/ , (4-v 2 -x)/ _ _3_4x m - l 8(x-x ) 



dx 2 x dx 



XT 



b 



where 



and 



9x 2 d s 



x = e y °/ 3 . 



For x < xq, 
For < x < x d , 



J = c 2 X- 2 +V / ^ + c 3 X- 2 - V^+*. 



J = C4 x- 2+ V^ + C5 x- 2 -^. 



(32) 

(33) 
(34) 
(35) 
(36) 



For Re(s) > 0, the boundary condition at y = —°° implies that C3 = 0. Integrating equation (32) 
from xq — £ to xq + £ yields the jump condition for df /dx at x = xo. After applying this jump 
condition and the continuity of / and r at x = Xd, we find that 



3x 



-v+\/v 2 +X -\ / v 2 +% 



where 
and 



b(hT - yx d U) 
T = I v (yx d )K v (y) - K v (yx d )I v (y), 



U = %(\\fx d )K v (\\f)-Kl,(\\fx d )I v (\\r). 
Since Iy(\\f)K v (\\f) - K^(\\f)I v (\\f) = 1 /\f, we find that 



(37) 

(38) 
(39) 



r(*=i) = - 



Since 



(/»> 



dll n T 



x=\ 



be 1 
~3~ 

ds 



(40) 



(41) 



x=l,s=0 



where n is a non-negative integer, we have 



(/»> = 



6 3"ci 



(42) 



5=0 
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Thus, 



1 + x f Xd 



4(v + l) 2v V-^o 
For fixed positive a in the limits — > 0, equation (43) gives 




(43) 



<<> = 2f>W (44) 

The full equation for (I 2 ) is very long and will not be quoted here. However, for fixed positive a in 
the limit xj — > 0, we find that 

(?) = 81 ^ + 6 ^ (45) 
K ' 4(3a + 2Z?) 2 (3a + 4Z?)' K J 

with Ls/h = {I 2 } /(I) given by 

Ls = 9(3a + 6fr) 

/ s 2(3a + 2Z?)(3a + 4^)' V 

From equation (43) (/) diverges for fixed xj as a — > (i.e. v — > 0). 

We now check the analytic results with a Monte Carlo solution of equation (20) with initial 
conditions P(y,l = 0) = 8(y — yo), and with yo = —10 and y d = —8. The Monte Carlo solution 
consists of iteratively incrementing a pair of numbers (l,y). During each step, we increase / by an 
amount 8/ = e 2y / 3 cp (or 8/ = e 2yd ^ 3 cp if y < y d ) and increase y by an amount 8v = acp ± k, where 
the ± sign is determined randomly with equal chance for either sign, and cp is a constant. As 
(p — > 0, each step becomes infinitesimal as is assumed in the Fokker-Planck equation. The value 
of k is chosen so that ((8j) 2 ) = 2£cp (i.e., k = ^/2b(p — a 2 (p 2 ). We stop incrementing / and y once 
y reaches 0, and record the value of / at the final step. We repeat this for 2000 ordered pairs to 
obtain (/) and {I 2 }. The results of this procedure with cp = 10~ 4 and a = 0.3 for various values 
of b are shown in figure 2, along with the analytic results of equations (44) and (45). We also 
plot Monte Carlo results with cp = 1, which provide a measure of the error associated with using 
a Fokker-Planck equation to describe the discrete stochastic process described by equations (16), 
(18), and (19). 

In figure 3 we plot Monte Carlo calculations of (/) as a function of initial field-line sepa- 
ration ro for b = 0.17, In/hi = 50, cp = 10~ 4 , and two values of a: a = 0.01 (solid triangles) 
and a = 0.29 (open triangles). The solid line is a plot of equation (44) for a = 0.01, and the dashed 
line is a plot of equation (44) for a = 0.29. The figure shows the increase in (Z) as a — > for 
fixed %d- 
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Fig. 2. — The solid line gives (/) from equation (44). The dashed line gives (I 2 ) from equation (45). 
The xs and squares give the values of (/) and (l 2 ) from Monte Carlo calculations with (p = 10~ 4 . 
The triangles and circles give the values of (/) and (l 2 ) from Monte Carlo calculations with (p = 1. 
For all data in plot, a = 0.3. 



5. Direct numerical simulations 



In this section, we study field-line separation in a static magnetic field using data obtained 
from direct numerical simulations of MHD turbulence. Galaxy clusters have very small mean 
magnetic fields, little rotation and thus little helicity, and very large magnetic Prandtl number P m , 
where P m = v/t), V is the viscosity, 7 and T] is the resistivity. Ideally, we would like to generate 
cluster-like magnetic fields self-consistently within a direct numerical simulation. However, at 
present it is not clear how to do this. Numerical simulations of turbulent dynamos in high-.P m 
plasmas driven by non-helical forcing find that amplified magnetic fields remain concentrated on 
very small (resistive) spatial scales (Maron & Cowley 2001). Yet the Faraday rotation produced 
by intracluster plasmas indicates that clusters possess considerable amounts of magnetic energy 
on large scales of order 1 — 10 kpc. We are thus faced with several alternatives. We can choose 
a numerical model that matches the helicity, mean magnetic field, and magnetic Prandtl number 
conditions of clusters, in which case the spatial scale of the model magnetic field is far too small. 
Or we can choose a numerical model that produces a large-scale field by, for example, including 



For clusters with anisotropic plasma viscosity, v is taken to be the parallel viscosity. 
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Fig. 3. — Monte Carlo calculations of (/) with b = 0.17, Is /Id = 50, cp = 10~ 4 and two values 
of a: a — 0.01 (solid triangles) and a = 0.29 (open triangles). The solid line is the analytic result 
[equation (43)] with a = 0.01, and the dashed line is the analytic result with a = 0.29. 

a mean magnetic field, using helical forcing, or choosing initial conditions that ensure that the 
magnetic field remains large-scale. 8 

In this paper, we follow the latter course. We initialize a simulation with a random-phase mag- 
netic field containing some amount of magnetic helicity, as well as random velocities with kinetic 
energies comparable to the magnetic energy. We allow the system to decay, leaving behind a large- 
scale magnetic field, which contains 10% of the maximum magnetic helicity for that magnetic 
energy. We then force the system non-helically at wavenumbers between 2n/Lb 0X and 4n/Lb 0X , 
where Z^ ox is the volume of the simulation cube. The forcing sustains a Kolmogorov-like spectrum 
of magnetic and kinetic energy. The result is a turbulent magnetic field that is dominated by large- 



8 Haugen et al (2003) found that non-helical dynamos with P m up to 30 result in a magnetic power spectrum Eb(k) 
proportional to k~ x . Such a spectrum is neither large-scale-dominated nor small-scale-dominated, since each loga- 
rithmic interval of order unity in fc-space contains the same amount of energy. However, they suggest that when both 
the ordinary Reynolds number R and P m are large, there is a fc -5 / 3 spectrum at large scales followed by a k^ 1 spec- 
trum at smaller scales, with the magnetic energy dominated by large-scale fluctuations. If this suggestion is correct, 
a very-high-resolution non-helical dynamo simulation with large R and large P m would be a self-consistent way to 
generate a large-scale magnetic field in cluster-like conditions. Such a simulation, however, would be beyond our 
current computational resources. 
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scale fluctuations, has a Kolmogorov-like inertial range extending to small scales, has zero mean, 
and has relatively little magnetic helicity, which we take to be a reasonable model for magnetic 
fields in clusters. We carry out a first simulation (simulation Al) on 256 3 grid points, and then 
use A 1 as an initial condiction for a higher resolution version (simulation A2) on 512 3 grid points 
with reduced resistivity and viscosity. We use incompressible simulations, which are a reasonable 
approximation for subsonic turbulence in clusters. The simulations are three-dimensional and pe- 
riodic, and the pseudo-spectral numerical method is described by Maron & Goldreich (2001). We 
use Newtonian viscosity and resistivity with P m = 1, since simulating P m 3> 1 in our isotropic- 
viscosity simulations requires a large viscosity that damps small-scale Alfven waves; small-scale 
Alfvenic turbulence is present in clusters due to anisotropic plasma viscosity (Goldreich & Sridhar 
1995, Quataert 1998), and plays an important role in field-line separation. In a companion paper, 
Maron, Chandran, & Blackman (2003), we study field-line separation in numerical simulations of 
different types of MHD turbulence and dynamo-generated fields. 

The simulation parameters for Al and A2 are summarized in table 1. In figure 4 we plot 
time averages of the magnetic power spectrum Et(k) [the total magnetic energy is / Et,(k)dk], the 
kinetic power spectrum E v (k) [the total kinetic energy is / E v (k)dk], and the total-energy spec- 
trum E tota [(k) = E v (k) + Eb(k) in simulation A2. We set 7t//g equal to the maximum of kEb(k), 
giving Ib = 0.25Lbox> and we set n/ld equal to the maximum of k 3 Eh(k). 



Simulation 


Grid points 


1(B) | 


Elm 


h/h 


a P m = v/ri 


Al 


256 3 





0.1 


23 


2.4 1 


A2 


512 3 





0.1 


50 


2.4 1 



Table 2: | (B) | is the strength of the mean magnetic field, H m is the magnetic helicity divided by 
the maximum possible magnetic helicity at that level of magnetic energy, Ib / Id is the ratio of outer 
scale to inner scale, a is the single-field-line diffusion coefficient in equation (5), and P m = v/t) is 
the magnetic Prandtl number, where V and T| are the viscosity and resistivity. 

We run each simulation until the power spectrum reaches a statistical steady state before we 
start analyzing field lines. For simulation Al, we take seventeen snapshots of the magnetic field 
separated in time by an interval 0AIb/u, where u is the rms turbulent velocity. To calculate Ls 
and (/) for an initial field-line separation ro, we introduce into each snapshot of the magnetic field 
2,000 pairs of field-line tracers whose initial separation vector ro is perpendicular to the local 
field. We use linear interpolation to obtain the magnetic field between grid points, employ second- 
order Runge-Kutta to integrate field-lines, and iteratively reduce the length step in the field-line 
integrations to achieve convergence. For simulation A2, we carry out the same procedure, but we 
use five snapshots of the magnetic field separated in time by an interval 0.2/g fu, and we use 20,000 
pairs of field-line tracers per snapshot. 
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Fig. 4. — Power spectra in simulation A2. The dotted line is the magnetic power spectrum Eb(k), 
the dashed line is the kinetic power spectrum E v (k), and the solid line is the total-energy spec- 
trum E tota i = E v (k) + E b (k) . 

We first seek to test the qualitative prediction of the Fokker-Planck model and previous theo- 
retical treatments (Jokipii 1973, Skilling et al 1974, Narayan & Medvedev 2001) that (/) and L s /l B 
asymptote to a constant of order a few as ro is decreased towards lj in the large-/g//</ limit. In 
figure 5, we plot (/) for simulation Al and simulation A2. The lower-resolution data of simula- 
tion Al suggest the scaling (/) oc ln(/£/ro) for 1^ < ro < 0.25/s, in contradiction to the theoretical 
treatments. On the other hand, for simulation A2, the curve through the data is concave down- 
ward for Id < ro < Ib- Moreover, figure 5 shows that the simulation Al data, and probably also 
the A2 data, have not converged to the high-Reynolds-number values of (/) for /g/16 < ro < Ib, 
values of ro that are within the inertial ranges (l^ to Ib) of both simulations. Also, the slope 
d{l)/d(\n(l B /ro)) for both ro < Id and ro < /s/10 decreases significantly when l B /ld is doubled. A 
comparison of the data for Al and A2 thus suggests that in the large-Zg/Z^ limit (/) asymptotes to a 
value of order several l B as ro is decreased towards l d , as in the Fokker-Planck model and previous 
studies. The same comments apply to the data for L$, which are plotted in figure 6. 

Because the definition of Is is not unique, we recalculate (/) and Ls, setting 27t//g equal to 
the maximum of kEt,(k), so that l B is twice its former value. Note that this affects both the unit for 
measuring distance along the field and also the distance to which field lines must separate. We plot 
the results in figures 7 and 8, which are qualitatively similar to the results based on our original 
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Fig. 5. — The average distance in units of l B that a field-line pair must be followed before separating 
by a distance l B , denoted (/), as a function of initial field-line separation ro for simulations Al 
(small circles) and A2 (large circles). 

definition of l B . 

We note that for ro = l B , (/) and L$ are by definition 0. The numerical- simulation data points 
that appear to be plotted above l B /ro = 1 actually correspond to ro just slightly smaller than l B , 
indicating that (/) and L$ are discontinuous at ro = l B in the numerical simulations. The reason is 
that for ro just slightly less than l B , some fraction of the field line pairs are initially converging and 
must be followed a significant distance before they start to diverge. 

We evaluate characteristic values of a and b, defined in equations (18) and (19), in simula- 
tion A2 by calculating the mean and mean-square increments to y = \n(r/l B ) for field-line pairs 
initially separated by a distance l B /S during a displacement of l B /4 along the magnetic field, using 
our original definition of l B [k/Ib = maximum of kEt,(k)]. We find that a = 0.29 and b = 0.17. 
These values are used to obtain the Fokker-Planck results plotted in figure 9. The Monte Carlo 
results in figure 9 use the same values of a and b and employ many small random steps (cp = 10~ 4 ), 
and are thus expected to reproduce the Fokker-Planck results. We also plot (/) for a random-phase 
version of A2, which is obtained by assigning each Fourier mode in simulation A2 a random 
phase without changing the modes' amplitudes. The figure shows that (/) is moderately larger 
in the direct numerical simulations than in both the random-phase data and the Fokker-Planck 
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Fig. 6. — Value of L$ as a function of initial field-line separation ro in simulations Al (small 
circles) and A2 (large circles). 

model. If we take a = 0.29, b = 0.17, lj = p,- = 43p e , and Is /Id ^> 1> the Fokker-Planck model 
yields Ls(Pe) — ^s(^d) — 4.5/g and the Monte Carlo model with order-unity random increments 
to \n{r/l B ) (i.e. 9=1) yields L s (p e ) ~ L s (l d ) ^ 6.5fe. 

In figure 10 we plot the probability distribution of the distance in units of l B that a pair of 
field lines in simulation A2 with initial separation ro = Id must be followed before the field lines 
separate by a distance l B . We define the function PDF(/) so that the probability that / lies in some 
interval is proportional to the corresponding area under the plotted curve. 

For clusters, Is /Id — te/Pi — 10 13 . In the \dxgQ-l B /l d limit, the numerical simulations and 
theoretical models indicate that L$ asymptotes to a value of order several l B as ro is decreased 
towards lj, and L$ is not expected to increase appreciably as ro is further decreased from lj = 
f>i to p e . Thus, Ls(pe) — Ls(ld)- To estimate Ls(ld) in clusters, we note that Ls(ld) — in 
simulation Al, and Ls(ld) — IO/5 in simulation A2. When the definition of l B is changed as in 
figure 8, so that 2n/l B corresponds to the maximum of kEb(k), then Ls(ld) — Hb in simulation Al 
and L$(ld) — 6.5l B in simulation A2. As mentioned previously, it is not clear which definition 
of l B leads to a more accurate prediction of Kj. We conclude from the direct numerical simulations 
suggest that Ls(p e ) — Ls(ld) — 5 — \0l B in the large-Zg/Zj limit. 
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Fig. 7. — The average distance in units of lg that a field-line pair must be followed before separating 
by a distance lg as a function of initial field-line separation ro for simulations Al and A2 (small 
and large solid circles). Same as figure 5, except that l B is redefined to be twice as large. 



6. Thermal conduction in time- varying turbulent magnetic fields 

In this section, we develop a phenomenology of particle diffusion in time- varying turbulent 
magnetic fields under the questionable assumption that the magnetic field is completely random- 
ized and reconnected on the eddy turnover time x at scale x = Ib/u, where u is the rms velocity 
and the velocity outer scale Iq is assumed equal to lg- A similar assumption was explored by 
Gruzinov (2002). For simplicity, we assume that the magnetic-field randomization occurs instan- 
taneously at regular time intervals of duration x. We assume x 3> A,/ V te> where v te is the electron 
thermal velocity and % is the Coulomb mean free path, so that particle motion along the field over 
a time interval x is diffusive. We assume that % <C Lg, so that particle motion along the field over a 
distance L$ is diffusive. We also assume that v te 3> u. 

There are three limiting cases. First, if x 3> L\/D\\, a particle escapes its initial field line 
through parallel motion and slow cross-field diffusion before the field is randomized. The "funda- 
mental random-walk step" is of length L$ along the field, as in section 2, and takes a time L\/D^. 
There are n steps = xDn/Lg such steps during each time interval x. The three dimensional distance 
travelled by an electron along the magnetic field during one "fundamental random-walk step" 
is ~ t/LsTb, and the three-dimensional distance travelled during n ste ps steps is ~ ^/n steps LslB- On 
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Fig. 8. — Value of L$ as a function of initial field-line separation r$ in simulations Al (small 
circles) and A2 (large circles). Same as figure 6, except that l B is redefined to be twice as large. 



the other hand, a fluid element travels a distance l B during a time x. Since L$ > l B and n ste ps 3> 1, 
the distance travelled by a single electron during a time x is much greater than the distance travelled 
by a fluid element, and the fluid motion can be ignored. The electron diffusion coefficient is then 
the same as in the static-field case. 

At the other extreme is the limit x <C l B /D\\. Since l B < Ls, a particle escapes a field line 
through field-line randomization before it escapes through parallel motion and slow cross-field 
diffusion, and the fundamental random-walk step is of duration x. The distance an electron travels 
along the field due to parallel diffusion, y/D\\%, is less than the distance l B that a fluid element trav- 
els. The net displacement of an electron (or ion) during one fundamental random step, denoted Ar, 
is then given by 

Ar ~ l B . {Al) 

Since the field is completely randomized after a time x, magnetic tension does not inhibit the 
wandering of fluid parcels over times x, as in Vainshtein & Rosner (1991) and Cattaneo (1994). 
Successive random steps are thus uncorrelated, giving a diffusion coefficient Ar^/x ~ l B /x = ul B , 
as in hydrodynamical turbulent diffusion. In this limit, the parallel diffusion of electrons plays no 
role. 

The third and intermediate case is l B /D» <^x<^ L$/D\\ . In this case, the field is again random- 
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Fig. 9. — The average distance in units of lg that a field-line pair must be followed before separating 
by a distance lg as a function of initial field-line separation ro for simulation A2 (solid circles), 
the random-phase version of A2 (open squares), Monte Carlo simulations with (p = 10~ 4 (solid 
triangles), and the analytic Fokker-Planck model (solid line). The Monte Carlo and Fokker-Planck 
solutions use a = 0.29, b = 0.17, and Ib/Ij = 50, values corresponding to simulation A2. 



ized before a particle can escape its initial field line through parallel motion and slow cross-field 
diffusion, and the fundamental random-walk step is of duration x. During one such step, a particle 
moves an rms distance ~ y^DyT along the field, which corresponds to an rms three dimensional 
displacement 

Ar~(l B ^D^y/ 2 . (48) 
The diffusion coefficient Ar 2 /x is then 

D ~ ^D\\ul B . (49) 

As an example, we consider the cluster A 1795 at a distance of 100 kpc from cluster center, 
with an electron density of 0.01 cm~ 3 and a temperature of 5 keV (Ettori et al 2002). We consider 
slightly superthermal electrons with Dy = Ks = 10 31 cm 2 /s. Churazov et al (2003) find evidence 
for turbulent velocities of order one-half the sound speed in the hot intracluster plasma of the 
Perseus cluster, which we take to be typical of A1795 as well, giving u ~ 350 km/s. We also 
assume lg = 10 kpc at this distance from cluster center, and take Ls = 61b- For these parameters, 
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Fig. 10. — Probability distribution of the distance (in units of /#) a field-line pair in simulation A2 
initially separated by a distance lj must be followed before separating by a distance /g. 

the plasma satisfies Z|/D|| < x < L|/D||. Equations (2), (13), and (49) then give 

K r ~ a/ k s w/ s , (50) 

which is about 0.3ks, or 3uIb- This value is roughly twice the estimate Ks/s/Ls, but given the un- 
certainties in both estimates, and the untested assumption that field lines are completely random- 
ized and reconnected on the time scale x, it is not clear that turbulent resistivity in fact enhances 
the thermal conductivity in clusters. 

7. Summary 

In this paper we consider the effects of field-line tangling and turbulent resistivity on the 
thermal conductivity in galaxy clusters. In the static-magnetic-field approximation, tangled 
field lines force electrons to move greater distances in traveling from hotter regions to colder re- 
gions, reducing Kt by a factor of ~ 5 — 10 relative to the Spitzer thermal conductivity Ks of a 
non-magnetized plasma for typical cluster parameters. It is possible that turbulent resistivity en- 
hances K T by a moderate amount relative to the static-field estimate for typical cluster conditions, 
but further work is needed to investigate this possibility. 
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